Decoherence and Entanglement Dynamics in Fluctuating Fields 
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We study pure phase damping of two qubits due to fluctuating flelds. As frequently employed, 
decoherence is thus described in terms of random unitary (RU) dynamics, that is, a convex mixture 
of unitary transformations. Based on a separation of the dynamics into an average Hamiltonian 
and a noise channel, we are able to analytically determine the evolution of both entanglement and 
purity. This enables us to characterize the dynamics in a concurrence-purity (CP) diagram: We 
find that RU phase-damping dynamics sets constraints on accessible regions in the CP plane. We 
show that initial state and dynamics contribute to final entanglement independently. 
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I. INTRODUCTION 

For newly emerging quantum technologies, robustness 
of quantum states is essential. In particular, stability of 
entanglement spread over a many-body quantum system 
is of fundamental importance for quantum cryptography 
and quantum information processing in general. A thor- 
ough understanding of processes known to destroy the 
desired quantum qualities — usually subsumed under the 
name of decoherence — is necessary. Phase damping rep- 
resents the case of pure decoherence, where coherences 
of a given state in a certain basis are subject to decay 
while probabilities, that is diagonal elements of the den- 
sity matrix, remain unchanged. It is well known that 
despite the rather simple nature of phase damping, it is 
clearly enough to disentangle quantum states [l], [2| ■ 

A widely encountered source of decoherence is growing 
entanglement between a quantum system and its envi- 
ronment. Yet, often decoherence is due to — or at least 
may be explained in terms of — random unitary (RU) dy- 
namics. Then, fluctuating classical fields are liable for 
the loss of quantum properties (sometimes also termed 
"random external fields" ) [1, 0] ■ Although it is known 
that RU dynamics is not the most general form of de- 
coherence d, [1] , it is of great practical importance. For 
instance, in quantum computers based on trapped ions, 
these classical fluctuations are believed to be the main 
source of decoherence Here, fluctuations are present 
both in the magnetic field of the trap and in the fre- 
quency of the laser addressing the qubits. RU dynamics 
moreover serves as a model to introduce decoherence in 
experiments in a controlled fashion Q. In the context 
of quantum error correction, RU processes are of special 
interest, because they represent the only type of errors 
that may be completely undone @. The present work 
will further reveal how decoherence of qubits due to RU 
dynamics allows for a surprisingly far-reaching analytical 
treatment. 

For a system of two qubits a possible means of charac- 
terization relies upon the relation between entanglement 
and entropy [l^iEl- The system's state q may be stud- 
ied in terms of a CP diagram, where entanglement (here 
measured through concurrence C) is plotted against pu- 



rity, P{q) — tr(f?^), a common measure for the mixedness 
of Q. This approach has been used in several theoreti- 
cal studies, where Bell states subject to single-sided dy- 
namics prj (by single-sided dynamics we denote the case 
where only one part of a bipartite system is affected), de- 
coherence in a random matrix environment 12], or two 
interacting atoms in a cavity [13| were considered. 

It was noted earlier that decoherence in ion trap quan- 
tum computers is not of simple exponential type 14], as 
would be expected from a Markovian master equation ap- 
proach. A theoretical model incorporating this fact thus 
calls for a more general treatment. From an axiomatic 
point of view (neglecting initial correlations) the dynam- 
ics of a quantum system is given in terms of completely 
positive maps (or quantum channels) In a Hilbert 
space of dimension d, these channels can always be writ- 
ten in terms of at most d"^ Kraus operators Ki such that 



(1) 



(throughout the article we denote the initial state by g 
and its map by g'). For RU channels there exists an ex- 
pression of the form ([T]) where every Kraus operator is 
proportional to some unitary operator, so that the dy- 
namics is a convex combination of unitary transforma- 
tions: g' = J2iPiUigUj, with Pi > and J^iPi = 1- 
Phase-damping channels stand out due to the require- 
ment to be diagonal in the distinguished basis, that is, 
they may be written in the form 



9mn — (fln|flm)(?n 



Drn.7i g^ 



(2) 



with__f|a„)} any set of d normalized complex vectors 
"system plus environment" models, the vec- 
tors I On) may be identified with environmental quantum 
states jlTf. For RU phase damping, a simple interpreta- 
tion of the |a„) is less obvious. 

In this article we study RU phase damping of two 
qubits. A fully analytical expression for both the en- 
tanglement evolution and the purity decay is obtained. 
Throughout the article the initial two-qubit state is as- 
sumed to be pure. The phase damping will be realized 
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as an ensemble average ({^U {t) gll^ {t)y^ , where the uni- 
tary (diagonal) time evolution may be regarded as arising 
from a stochastic Hamiltonian of the form 

H{t) ^ oji{t)a^ t + uj2{t)t + uj3{t)<T^ (S) . (3) 

With (T^ we denote the diagonal Pauli spin operators for 
qubits A and B, respectively, the time dependence is due 
to the stochastic processes uji{t),uj2{t), and u}3{t). In 
an experimental realization using two-state atoms this 
would correspond, for example, to a fluctuation of the 
individual Zeeman levels described by iOi{t) and ^2(^)1 to- 
gether with an instability in the interaction of the atom's 
energy eigenstates, described by (^^{t). Such a Hamilto- 
nian is used to describe the spin-spin interaction in nu- 
clear magnetic resonance (NMR) systems 1^] . Also note 
that the third term of Eq. ([3]) describing the interaction — 
later referred to as pure two-qubit phase damping — is im- 
plemented in the realization of a phase gate in recent ion 
trap experiments Q- 

Recently, for a quantum system subject to single-sided 
dynamics a very simple evolution equation for entangle- 
ment was found [19]. This evolution equation allows for 
a factorization into a functional characterizing the quan- 
tum channel and the amount of the initial pure-state en- 
tanglement. In our case here, the evolution will in general 
be nonlocal. Yet, we are also able to give some exact evo- 
lution equations for entanglement under phase damping, 
where the channel and the initial state enter indepen- 
dently. 

The article is organized as follows: In Sec. II we dis- 
cuss our model of RU phase damping, where we ben- 
efit from a separation of the channel into a reversible 
part due to a mean Hamiltonian and an irreversible noise 
channel. This separation enables us to analytically study 
changes in entanglement and purity. In Sec. HI we dis- 
cuss RU single-qubit channels in the context of the CP 
diagram. Here we are able to give an alternative ex- 
planation of the results obtained in Ref. [ll| , where for a 
maximally entangled state subject to unital (mapping the 
completely mixed state onto itself) single-sided channels 
certain bounds to the accessible area within the CP dia- 
gram were found. Our approach, which employs the so- 
called Jamiolkowski isomorphism, will also be well suited 
for explaining part of the findings at hand. In Sees. IV 
and V we analyze in detail the evolution of entanglement 
and purity of a pure separable and a general pure initial 
state under RU phase damping, respectively. 



II. RU PHASE DAMPING 

The phase damping we consider shall be realized as an 
ensemble average g' — ({^U (t) gll'' (t))) , where we may set 
U{t) = e~^-fo d.TH{T)^ ^ unitary map based on a stochas- 
tic diagonal Hamiltonian H{t). As an example let us 
consider a single-qubit channel, where a generic diagonal 
Hamiltonian may simply be set to H{t) = uj{t)<Tz with 



some stochastic process uj{t). Assuming ijj{t) to add up 
to a total perturbation of the qubit, the central limit 
theorem tells us that VI :— J* dTUj{T) is a Gaussian ran- 
dom variable that is determined by its mean, /i := 
and its variance, cr^ =: {{^'^)) — M^, hence leading to 

Substituting fi 6 and i(l — e^^"' ) — )■ p, we thus 
arrive at the most general form of a single-qubit phase- 
damping channel 

^ - [{{e^^'))g2i g22 ) 

= e-'i^^{{l~p) g + pa,ga,)e'i''^. 

Thus, single-qubit phase damping may always be written 
in terms of RU dynamics [5] . From our example we see 
furthermore that it suffices to consider Gaussian fields. 

For two qubits the corresponding stochastic diagonal 
Hamiltonian is 

H{t) = uji{t)a^ ® t + uj2{t)t ® + uj3{t)af ® af . 

Irrespective of possible correlations among the accumu- 
lated phases (iTWfc(r) =: fi^, we find that the phase- 
damping channel can be decomposed according to 

g' - U^{D*g)Ul (4) 
= D^iU^gUl) 

into a unitary, reversible part [/^ based on the mean 
Hamiltonian iJ^ = ^laf 1 + (^2^^ cr^ + ^J■3C^^ crf 
ifJ-k ■= ii^k))) and a nonunitary, irreversible noise chan- 
nel D (see Appendix El for more details). For brevity of 
notation we here made use of the Hadamard poduct ★ 
of matrices, that is, the pointwise multiplication of two 
matrices of the same size [l^l . 

In close analogy to the single-qubit case we can further 
deploy the RU phase-damping channel by assuming n to 
be a Gaussian process. It should be noted, however, that 
our analysis could easily be extended to more general 
statistics. We can then use the characteristic function 

m 

where fl :~ ((fi)) denotes the mean value, S is the co- 
variance matrix, and Efe/ := (((fifc — i-j,k){^i — ^J■l))) for a 
Gaussian process fl. 

With no correlations between the stochastic processes 
uji{t),uj2{t), and W3(t), the covariance matrix is diagonal: 
S = diag(2<;J, 2c:|, 2^1). Note that -\/2c:fc simply denotes 
the standard deviation of the Gaussian random variable 
flk ■ The preceding commutativity relation @ now allows 
for a separation of the channel into its mixing dynamics 
and its entangling dynamics: Any change with respect to 
the purity of the system is due to the noise channel D, 
while any increase in entanglement can only be evoked 
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by (the interaction part of) the mean Hamihonian H^. 
Of course, a decrease in entanglement might as well be 
due to D. However, any gain in entanglement may only 
be attributed to the unitary map C/^. This separation 
enables us to calculate changes in purity and concurrence 
separately. 

In order to describe the loss of purity of the system, 
we need to account for the action of the noise channel, 
b, only, 

D-kQ ^ PiQ + P2{cr^ ® Q ® af) (5) 
+ p^{(Tf ®t gaf® 1) +P4(1 gt® erf), 

where the probabilities pi with '^pi — 1 are determined 
by the variances of the stochastic processes (see Ap- 
pendix |B|. The purity for the final state is easily ob- 
tained and equates to 

P{b*g) = {p\+pI+pI+pI) (6) 

+ 2(piP3+P2P4) |(fT^)|' 
+ 2(piP4+P2P3) |(af)P 

+ 2(pip2+P3P4) |(a^(8af)|2. 

Here and in the following we write ( • ) for expectation 
values with respect to the initial state |V'o)- While the 
equation for purity decay under RU phase damping is 
valid irrespective of further details of the initial state, 
the change in entanglement will be more sensitive to the 
initial conditions of the two-qubit system. In order to 
calculate the entanglement created or destroyed in the 
phase-damping process, we thus need to make some fur- 
ther assumptions about the initial state. 

III. RU SINGLE-QUBIT CHANNELS IN THE 
CP DIAGRAM 

In the analysis of RU phase damping we want to fol- 
low the example of Ref. where single-sided unital 
dynamics was analyzed with respect to entanglement de- 
cay as a function of purity decay. The entanglement for 
a two-qubit system can be easily computed in terms of 
concurrence, C{g) = max{0,Ai — A2 — A3 — A4}, where 
Af > A2 > A3 > \\ are the eigenvalues of the positive 
matrix R = g {oy ® ay) g* {ay (8) ay) [22]. The purity P 
is defined as the trace over the squared density operator: 
P{g) = ti{g'^). Accessible values of purity P and con- 
currence C are not independent of each other. Rather, 
it has been shown that — depending on the entanglement 
measure used — there exist states that maximize the re- 
spective entanglement measure at a given mixedness, also 
called maximally entangled mixed states (MEMS, see (TT| 
and references therein). This fact is nicely illustrated in 
the concurrence-vs-purity diagram (CP diagram; see Fig. 

Aside from the border of all physical states given by 
Cmems; the diagram reveals the accessible region (blue) 
for initial Bell states under single-sided unital channels 
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FIG. 1: (Color online) The CP diagram. The hght gray area 
shows the region of all physical states of two qubits, bounded 
by the maximally entangled mixed states (Cmems). The pos- 
sible C{P) relations for a Bell state subject to single-sided 
unital dynamics or, equivalently, of all doubly chaotic states, 
is highlighted in blue (dark gray). This region is bounded 
by the curves valid for Werner states Cw and for Bell states 
under single-sided phase damping Co, respectively. The dia- 
gram also reveals the maxium C{P) relation of doubly chaotic 
states with Bell rank 3, Cs{P) (dashed line). 



f (X) 1, that is, unital channels acting on either one of 
the subsystems, only. This accessible region found is 
bounded by the curves Cw and Co that are realized for 
single-sided depolarization and single-sided phase damp- 
ing, respectively. The former is also valid for the so-called 
Werner states gw = P \4's) {4's \ + (1 — p)l/4, a con- 
vex mixture of the pure and maximally entangled singlet 
state \tps) = (|01) — |10))/V2 and the maximally disor- 
dered two-qubit state 1/4. 

These borders will also play an important role in our 
findings. We start by rederiving the results found in 
Rcf. 11] in an alternative way. It is known that in case 
of single-qubit dynamics there is a one-to-one correspon- 
dence between unital channels and channels that are RU 
0. By applying the quantum channel £ to one part 
of a maximally entangled bipartite state, the so-called 
Jamiolkowski isomorphism 24] introduces a duality be- 
tween quantum channels defined on a Hilbert space H 
and quantum states living on a Hilbert space % ® H 
of squared dimension. It follows straight away that the 
Jamiolkowski state of a RU channel is doubly chaotic [Ij , 
that is, the state has maximally disordered subsystems: 
trig = tr2£i = l/d (here, tr^, i = 1,2, refers to a par- 
tial trace over the first and second subsystems, respec- 
tively, while d — dimH). In this vein, any single-qubit 
RU channel may be identified with a doubly chaotic two- 
qubit state. Any such state may in turn be obtained by 
applying a local unitary transformation [/ U to a Bell- 
diagonal state [15, that is, a state diagonal in the basis 
of BeU states 
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with J2iPt = ^- 

For such a state the accessible area in the CP diagram 
may be determined easily (see Appendix[C]). As a general 
lower bound we find the relation C{P) = \/2P — 1 — Cd, 
whereas the upper bound depends on the Kraus rank k 
of the corresponding qubit channel £, that is, on the 
minimum number of Kraus operators needed in ([1]) . For 
fc = 2 we have 



C2iP) = C^ = V2P~^, P > 



1 

2' 



(7) 



again identical to the lower bound in the CP diagram, 
whereas for fc — 4 wc have 

C4{P) = C^^\ (Vl2P-3-l) , ^>^, (8) 

reproducing the upper bound. In addition, we get the 
relation 



C3(P) = i(V24P-8-l), P>| 



(9) 



giving an upper bound of the allowed area for single-sided 
RU channels of Kraus rank 3 (see Fig. [1] dashed line) . 

Thus, we arrive at the same results that were obtained 
by applying single-sided unital channels to a pure and 
maximally entangled two-qubit state We find that 
the accessible region in the CP plane depends on the 
Kraus rank of the RU channel. It is then quite obvious 
that the lower and upper border are given by single-sided 
phase damping (Kraus rank 2) and depolarizing (Kraus 
rank 4) , respectively. The various bounds in the CP dia- 
gram may be attributed to a characteristic trait of doubly 
chaotic states. This is also true for the behavior discov- 
ered in where for an initial Bell state of two non- 
interacting qubits with symmetric coupling to a cavity 
with photon number n — ^ oo a convergence to C3 (P) was 
found. In this limit, the Bell state is simply mapped to 
a doubly chaotic state with Bell rank 3, hence explaining 
the asymptotic behavior. 

For random local unitary (RLU) channels, that is, 
channels of the form £{q) = '}2,r,sP-rs{Ur®Vs)Q{Ul®V}), 
with g Prs = 1 and unitary U^^Vs, it is easy to see that 
a Bell state is mapped to a doubly chaotic state. Com- 
parison with Eq. ([5]) shows that the action of the noise 
channel D is exactly of this RLU nature. Therefore, it is 
by no means surprising that the bounds play a role in our 
findings. Note, however, that our studies cover not only 
the case of maximally entangled initial states, but of pure 
initial states in general. Based on our derivation it also 
follows immediately that the accessible region for bilocal 
unital dynamics Si ig) £2 — a question raised in Ref. — 
is in fact identical to the area obtained for single-sided 
unital channels. 



IV. SEPARABLE INITIAL STATE 

In Sec. [n] we discussed the possible separation of 
the RU phase-damping channel, making it possible to 



separately calculate the evolution of purity and entan- 
glement. The general formula for purity of the final 
state was already given, Eq. (j6|), whereas in order to 
study entanglement, we need to make further assump- 
tions about the initial state. In this first section we 
take the initial state to be separable. Due to the com- 
mutativity of the noise channel D and the unitary part 
L/^, we can start out by bringing a pure separable state 
g := \ipo) (V'ol := I'/'i) (0i| ^ |02) (02| into the mixture 
D-kg, prior to applying ?7^. For the final state g' we will 
then calculate the concurrence C. For generic single- 
qubit states we can let 



, z-1,2. 



with ai,f3i e C,|aip -I- = 1. The calculation of 
the concurrence of g' involves the determination of the 
eigenvalues of the positive, non-Hermitian matrix R = 
g'{<jy (g) ay)g'*{ay ® ay). We find (see Appendix ID]) 



C{g') = sdV-o)) • f{£), 



(10) 



4|ai||/3i||a2||/32| - y^var(a,^)y^var(af), 
niax|o, {pi - P2) sin ^3 

(P3 + Pi)"^ sin^ ^3 -I- 4p3P4 cos2 ^aj. 



we denote the variances of 



where 
5(1^0)) 

fin 

With var(cr*) 

a\ {i = A,B), which can be taken with respect to the 
initial state IfAo)- The mean value /i3 = (($^3)), as well as 
the probabilities pi, were introduced in SecHIland refiect 
mean value and variance of the fluctuating fields (see also 
App.ll). 

The concurrence thus allows for an evolution equa- 
tion factorizing into two functionals / and g that depend 
on the phase-damping channel £ and the initial state of 
the qubits IV'o): respectively. The evolution equation ac- 
counts for the gain in entanglement due to the unitary 
map C/^ through the function /, while g simply gives a 
scaling factor for the amount of entanglement a certain 
initial state may achieve. 

When trying to detect the accessible area within the 
CP diagram we need to know the maximum of the con- 
currence with respect to the entangling part of the chan- 
nel. This maximum is given for sin //3 = 1 or, equiva- 
lently, ^13, G {{2k + | k G Z}, which leads to the 
simple equation 



= max 



{o,.g(|^o)) {pi 
{0,5(1^0)) 



P2 - P3 - P4 



)} 



(11) 



-{^1+4) 



- 1 
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A. Pure two-qubit phase damping 

We want to discuss the results obtained so far in a set- 
ting where we admit only pure two-qubit phase damp- 
ing. By this we mean the situation where the phase 
damping is caused by a fluctuation in the part of the 
Hamiltonian describing the coupling of the two qubits 
only. Note that this setting is also of relevance in the re- 
alization of a phase gate in recent ion-trap experiments 
The Hamiltonian then consists of the single term 
H{t) = u!3{t)a^ €5 erf . Recall the definition of the vari- 
ances 2<;^ = {{nl)) ~ {{nk)f (cf. Sec. mi). We may there- 
fore let the variances <?i = ^2 = 0, and Eq. (fTT]) simplifies 
to 



whereas the purity simplifies to 
1 



P{9') 



:(1 
1 



(12) 



(13) 



(l-e-'^^)-\{af)\'-\{a 



The simple form of both the maximal concurrence and 
the purity enables us to give the former as a function of 
the latter: 



tA\2\ 



1 



rB\2) 



'2P - 1 



TA\2lrrB\2 



(14) 



we get CMax(P) 



N ote tha t for (cr^) = {af) 
\/2P — 1, which is exactly the equation obtained for the 
lower bound Co shown in Fig. [TJ 

In Fig. [2] we show the CP diagram where we high- 
light the excluded area. In order to show exemplary dy- 
namics in the diagram we let both mean and variances 
of the stochastic processes in Eq. ([3]) depend quadrat- 
ically on some "time parameter" t. This introduces a 
one-parameter class of channels 8t mimicking continuous 
dynamics. A possible evolution is shown as blue line. 

Equation (jl4[) relating maximal concurrence and pu- 
rity allows for a representation of the accessible region 
under pure two-qubit phase damping (for an initially 
separable state). For this we define Z as the expec- 
tation value of the number of excited states, that is, 

Z := 1 + - , ' . When maximizing the function 
C-MB.,i{P) for fixed P and Z, we get the requirement 
(a1^) = (a?) =: (a,). Inserted in (Illl), this leads to 



Cm..(P,(ct,)) 




(^.>^)- (15) 



This relation can now be visualized in the concurrence- 
purity- (cTz) space [cf. Fig.[2](b)]. This representation was 
also used in |26j . where the physically allowed region of 
a two-qubit system with respect to concurrence, purity, 
and energy was studied. 
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FIG. 2: (Color online) (a) For a pure separable initial state 
subject to pure two-qubit phase damping, the maximally 
achievable concurrence is bounded from above by Cn- The 
corresponding forbidden area is highlighted in red (dark gray) . 
The blue (dashed) line represents the evolution under a one- 
parameter class of phase-damping channels (see text for de- 
tails), (b) The accessible value of concurrence depends on 
expectation value {ai'^) evaluated for either one of the qubits 
i = A,B. Only for (cr^) = (cr^) = may the relation Cu 
(red, thick line) be obtained. 



B. General uncorrelated phase damping 

Next we want to consider the general case of uncorre- 
lated phase damping. Remember that by uncorrelated 
we mean a diagonal covariance matrix E. Here, we use 
the full dynamics, that is, the part of the Hamiltonian 
describing pure two-qubit dephasing as well as the parts 
acting locally on both qubits. 

As in the preceding section, the concurrence attains its 
maximum for (cr^) = (ct^) = 0, leading to the identities 



C 



P 



Note that for identical variances <;i — — 'is —'■ this 
implies C = i (ae-^":" - l) and P = i (Se^'*'^' -I- 1) , which 
is easily transformed into the relation 

C(P) = i(Vl2P-3-l), P>i (18) 
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FIG. 3: (Color online) (a) In case of general uncorrelated 
phase damping, the maximally achievable concurrence is iden- 
tical to the value given by the Werner states, Cw As before, 
the forbidden area is highlighted in red (dark gray). The role 
of correlations may be observed in (b), where both bounds 
are violated. In both diagrams, an exemplary evolution un- 
der a one-parameter class of phase-damping channels is shown 
(dashed blue lines). 

which coincides with the relation obtained for the Werner 
states, . 

Furthermore, from Eqs. and ([T7)) we conclude 

that C'^ + C + 1 < 3P, so that the case of all vari- 
ances being equal already yields the maximal relation of 
concurrence as a function of purity. We conclude that 
Cmb.x{P) = (-\/12P — 3 — l)/2 in case of uncorrelated 
phase damping. For pure, separable states subject to 
RU phase damping, we can therefore identify a forbid- 
den area in the CP plane, the border of which is given 
by the Werner states Cw In Fig. [3] (a) we highlight the 
excluded area for general, uncorrelated phase damping. 

Note that these results are only true for uncorrelated 
phase damping. In Fig. [3] (b) we show an example of 
correlated phase damping, that is, with nondiagonal co- 
variance matrix E (see Sec. |lll. This accentuates the 
fact that the borders in the CP diagram are valid only 
in the uncorrelated case. We ascribe this violation of the 
bounds in case of correlated phase damping to the fact 
that the irreversible part of the channel D is in general 



no longer RLU (see Sec. IIIII) . 

V. GENERAL PURE INITIAL STATE 

In this section we will analyze the evolution of a general 
pure initial state under phase damping, that is, we drop 
the restriction of separability for the initial state. For our 
purposes, that is, for entanglement and purity evolution 
under phase damping, the most general (pure) intial state 
reads 

leaving only four parameters (a, x, ^2) to fully charac- 
terize the initial state of the two-qubit system (see Ap- 
pendix |E]) . 

While Eq. ([6]), giving the purity of the final state, is 
still valid, we need to estimate the decay of the initial 
concurrence when the state is subject to the mixing part 
of the phase damping channel, D. In close analogy to 
the method described in Appendix |D1 we are again able 
to estimate the concurrence analytically: 

C{g) = Max|o, 

\/ (P? + 2piP2kP - 2pip2\s^ - r2| - 

{pI + pI)\s\^ + 2p3P4|r|2 + 2p3Pi\s^ - r2|}, 

where we have defined s := {ipol T,y and r := 
(■0o| EzEj; \iPq). Here we have also used the abbreviations 
Ej := CTi (g)cri {i = x,y,z). 

Together with the purity (|6]) we can now examine the 
possible C{P) combinations in more detail. In order to 
get illustrative results, we again have to put some further 
constraints on the dynamics. 

A. Single-sided phase damping 

For a single-sided phase-damping channel, that is, act- 
ing on either one of the qubits, we can generalize the re- 
sults obtained for maximally entangled states in Ref. [ll| 
to the case of pure initial states with arbitrary entangle- 
ment. Note that for this scenario the channel is simply 
given by £>' = (!- q)g + q{a^ l)e{o-f «) 1). 

For single-sided phase damping it is easy to arrive at 
the evolution equation of entanglement, which is given 
by 

C{g') = f{E) ■ Co, 

where f{£) = 1 — 2q. The state-dependent behavior 
within the CP diagram thus stems from the evolution 
of purity, which is given by 

Pig') = l-2q{l-q) [l-\{af)\^], 



0.2 0.4 0.6 0.8 1.0 

purity 

FIG. 4: (Color online) Accessible region for a single-sided 
phase damping channel for states of varying initial entangle- 
ment Co = 0.2 (0.4, 0.6, 0.8, 1). For Co = 1, the lower and 
upper bounds are equal, given by the known curve Co. 

where |(cr^)|2 = cos^ 6*1(20 - 1)^. 

We find that the C{P) relation for given initial entan- 
glement now only depends on the angle 9i . The extremal 
cases may be identified as 

r(P\ = l Vgg(2^-1) ifcos^gi^O 
^ ' 1 ^Cl + 2{P-l) if cos2 Oi = 1 

where Cq denotes the concurrence of the initial state 
I'f/'o)- For all other cases (0 < cos^ 6'i < 1), one can see 
that the relation C{P) is bounded by these two curves. 
For maximally entangled initial states the two functions 
coincide, again giving the by-now well-known relation 
C{P) = V2P- 1 = Cd [ni. For nonextremal initial 
entanglement < Co < 1, however, the two curves part, 
bounding the accessible region from above and below. 
For states of arbitrary initial entanglement subject to 
single-sided phase damping we conclude that there exists 
a maximal value of purity, P< = 1 — Co/2, above which 
all states are still entangled. Again, the results are easily 
illustrated using the CP diagram (see Fig. 0]). 

B. Pure two-qubit phase damping 

In case of pure two-qubit phase damping, the noise 
channel is of the simple form £»' = (! — q)Q + 9(0"^ ® 
af)g{a^ ® af). Again we are able to give the evolution 
equation of entanglement, 

C\g') = Cl-2q{l-q)[\r\^~\s\'~\r^-s^\\ 

Cl~f{£)-g{\M). (19) 

where now f{£) = 2q{\ — q) and g{\iljQ)) = |rp — |sp — 
— I . Entanglement depends on the initial state of the 
two qubits through the functional g. Note that here the 
evolution equation no longer follows the factorization law, 
as in the case of single-sided channels. Yet the evolution 
is fully determined through the functionals / and g. The 



0.2 0.4 0.6 0.8 1.0 

purity 

FIG. 5: (Color online) For pure two-qubit phase damping 
we can identify states with robust entanglement. The plot 
is obtained for an initial state parameterized by x = 0, Si = 
37r/16, 62 ~ vr/8 (see text), leading to robust entanglement for 
Co < C> « 0.5142 (dashed red line). The blue (thick) lines 
show exemplary evolution of pure states with initial concur- 
rence Co = 1.0 (0.8, 0.6, 0.4, 0.2). 

purity equates to 

P{g') = 1 - 2g(l -q)[l- \{a^ ® af )p] (20) 
with |((T^ (g) erf )p ~ [cos 01 cos 6*2 — sin^i sin02 • cosx • 

Eq. ([TO)) highlights the possibility of robust entangle- 
ment under pure two-qubit phase damping. It is easy 
to see that for g{\tpQ)) — concurrence remains con- 
stant, whereas purity [Eq. ([20]) ] may, in general, still 
decline. We find that this is the case for states with 
X € {fc7r|fc e Z}, only. The square of the concurrence 
then becomes 

, , _( s^- 49(1 - q)[Ci - r2] if s2 > 

The necessary condition for robust entanglement under 
pure two-qubit phase damping is thus given by < 
and X — ki:. To be specific, the necessary conditions for 
X — k2TT [x = {2k + 1)tt] are given by 

g2 ^ j,2 

<^ 

Co < C> :^max(h] ""^^ "''^^^li , [-] 

^ — ^ I ^1 P2 + 1 COS ^1 COS P2 — 1 I 

hereby defining the upper bound C>, such that states 
with initial concurrence Co < C> exhibit robust entan- 
glement. Note that for cos 6*1 cos 6*2 — ±1, s'^ — and, 
thus, constant concurrence C follows immediately. In 
Fig. [5] we show an example of robust entanglement. It can 
be clearly observed that for initial concurrence smaller 
than C> and despite the loss of purity, the entanglement 
is persistent. Beyond the observation of robust entangle- 
ment we also see that, while the map of the maximally 
entangled state follows the line Co , as would be predicted 
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(the Kraus rank of pure two-qubit phase damping is 2, 
see Sec. Ill), for states with initial concurrence Cq < 1, 
the restricted area no longer plays a role. 

VI. CONCLUSIONS 

We have studied a two-qubit quantum system subject 
to phase damping caused by fluctuating fields. Irrespec- 
tive of possible correlations between the individual con- 
stituents of the introduced stochastic Hamiltonian, the 
quantum channel was shown to allow a separation into 
two parts, representing reversible and irreversible parts 
of the dynamics, respectively. The separation enabled 
us to estimate the evolution of the system's entangle- 
ment, as well as purity, analytically. When combined, 
the results allow us to identify exclusive regions within 
the concurrence-purity plane, suggesting a connection be- 
tween RLU dynamics on a Bell state and phase-damping 
processes on a separable initial state. These forbidden 
places, however, do not play a role in the evolution of 
pure states of arbitrary initial entanglement < Co < 1. 

Our approach enables us to generalize results obtained 
in Ref. [ll| , where the action of single-sided unital chan- 
nels on Bell states was studied. Here we are able to 
analyze the effect of single-qubit phase damping on pure 
two-qubit states of arbitrary initial entanglement. For a 
certain class of phase damping we have identified neces- 
sary conditions leading to a robustness of the entangle- 
ment present in the two-qubit system. This robustness 
is closely related to the entanglement of the initial state; 
more precisely, depending on the initial state, we can give 
an upper bound C> such that for all states with concur- 
rence Co < C> the entanglement is robust. The crucial 
conclusion to be drawn from these findings is that in or- 
der to preserve entanglement, sometimes less is more. 
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Appendix A 

In Sec. In] we present the decomposition ^ of the 
phase-damping channel into unitary part and noise chan- 
nel. In order to see this, recall the definition of the RU 
channel 

g' = {{U{t)gUHt)))^ ^ 
The unitary map is given by U(t) — e ^•'o '^'^^i'^) ^ where 
in our model 

H{t) = LOi{t)af «) 1 + L02{t)t ® erf + W3(t)cr^ ® crf . 

A simple calculation gives the diagonal map 

o' - //e'^'=i Wo 

gran ~^ \\ ^ // fmn 

— • ^rrin^ningniTLi (-^-^) 

where we have used the compact notation c™" := (e,„ — 
e„)fe with vectors ei = (0,0, 0), 6*2 — (0,1,1), 63 = 
(1, 0, 1), 64 = (1, 1, 0), as well as the abbreviations = 

e'ELi--{(f^.» and Ann = ^^e^^ZLi J)^ 

The diagonal map Z^^Jm may be rewritten in terms of a 
diagonal unitary transformation f/^ = e"*^*"* such that 
= iUtj,),nmiUtj,)nm whcrc the Hamiltoniau iJ^ is de- 
termined by the mean values {{^k)) l^k '■ 

Hf, fliaf (g) 1 + ^2! «) crf + fi3<7f «> crf. 

Commutativity of unitary part and noise channel D 
is then a direct consequence of the diagonality (|Aip of 
the two maps. Formally, we may thus write 

g' - U^iD^g)Ul 
= D^'-kibi.g) 
= Di.iDf'^g) 
= b^{U,gUl). 

Appendix B 

The probabilities used in Sec. |lT] in terms of the vari- 
ances 2i,1 — {{^k)) ~ ((^fe))^ given as follows: 
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Appendix C 

In Sec, mil we argue that any unital single-qubit chan- 
nel may be identified with a doubly chaotic two-qubit 
state. A given doubly chaotic state of two qubits may in 
turn be obtained by applying a local unitary transforma- 
tion [/ (g) y to a Bell-diagonal state [l^ : 

with probabilities pi, {|5'f )}i=i,...,4 denotes the basis 
of Bell states. By rearrangement of the Bell basis let 
Pi > P2 > P3 > P4, so that m denotes the Bell rank of 
g, that is, the minimum number of Bell states needed to 
represent the state g. Note that the Bell rank is equal 
to the so-called Kraus rank of the channel, which gives 
the minimum number of terms in the Kraus represen- 
tation ([T]). Using the Bell-diagonal representation, both 
concurrence and purity are of a very simple form: 

m 

=niax{0,2pi-l}, P{q) =Y.pl (CI) 

4=1 

The maximum relation C{P) now depends on the Bell 



rank m. It may be obtained by minimizing the purity for 
given concurrence. Let pi —: 1 — then from Eqs. (|Cip 
we can conclude that the purity is minimal if the remain- 
ing weights Pi with i > 1 are equal, pi = qjim — 1). We 
thus obtain 

C{g) = max{0, l-2g}, 

/ Tfl \ 

Pig) > l-2q[l-q , 

and we can immediately give an upper bound for the con- 
currence as a function of purity: C{P) < 2 ^^^^^ \l' ^^m-l ~ 

2^^^^ =: Cm{P). Inserting the nontrivial Bell ranks 
m = 2,3,4, we get the upper bounds of the accessible 
C{P) relations for single-sided RU channels with corre- 
sponding Kraus ranks. Moreover, we see that for ar- 
bitrary Kraus rank fc, C{P) is bounded from below by 

C2{P). 

Appendix D 

The calculation of the concurrence involves the de- 
termination of the eigenvalues of the positive, non- 
Hermitian matrix 



n = g g 

e""^ IV'o) I - jp?sin/i3 (^/-ol S^e"*"^ -piP2COS^3 (i/'o| S^E^e""^! 
+e-"^S, { -P1P2COSM3 {%\^ye-"^ - ip^sin/ig (V-ol S,,E,e-''^} 
+e-'^a, ® 1 IV-o) {ip3sin//3 Y.y{a, ® t)e-"^ +p3PiCOSH3 (V'SI (g) cr,)e-"^^ 
+e-"^l (g) a, IV-o) {p3P4COSAi3 (V'SI ^yi'^z «) 1)6'"^ + iplsin^is (V'ol (g cr.)e"*'^} , 

I 

where we have defined :— Ui ® Ui {i = x^y,z) and We find the existence of two i?-invariant lin- 
(p := lij,T,z/2. ear subspaces S, S C C, spanned by the vectors 

{e-"^|V'o),e-"^S.|^o)} {^i.^^t} and {e-"^a, (g 
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Iji/io), e"*"^! (X) (TzIV'o)} =: {1^1,1^2}, respectively. For the 
determination of the eigenvalues of R we thus have to 
find a basis of orthogonal eigenvectors of spaces S and 
S. For instance, for eigenvalues A^, in S, we need to 
consider 

R{i'i + ki>2) — 0,1^1 + i>V2 + k {cvi + di'2) 
= A^ {vi + kv2) 

where a, b, c,d,k € C, leading to the quadratic equations 
A| = ^ ± ^ + be. Note that in the formula 

for the concurrence there is a sum of the square roots 
of the eigenvalues of the matrix R. Therefore the rela- 
tion (A+ ± A_)^ — (a + d) ± 2\/ad — be will prove to be 
quite useful. For the second subspace S, the situation 
is analogous and the eigenvalues of R are thus given by 
the set {A^, A^, A^, A^}. The concurrence is eventually 
given by 

C{g') = Max{0,4|ai||/3i||a2||/32| (a+- A_ - A+-A_)} 

= g{\4'o)) ■ f{£), 
where we define 

/(£) Max|o, (a+ - A_ - A+ - A_) | 
= Max|o, {pi - P2) sin ^3 - 

\J {p?.+Piy sin^ + 4p3P4 cos2 /i3 1 



and 

g(IVo)) := 4|ai||/3i||a2||/32|. 



Appendix E 

By use of the Schmidt decomposition, any pure two- 
qubit state may be written in the form j?/') — Ui (S) 
(\A|00) + yr^|ll)), where UuU2 G SU{2) (j. 
Note that the concurrence depends on the single pa- 
rameter a only and equates to C{\ip)) = 2^a{l — a). 
Using the Euler angles {ipi,Oi,Xi) the unitary rotations 
may be written in the form Ui = e~i'^^'^'e~i^''^''e~i^''^' , 
i ~ 1,2. Due to diagonality of both the phase-damping 
channel and the last rotation, e~2'^i<^^ e~'^'^^"'% we 
may reverse their order (diagonal operators commute). 
When interested in entanglement and purity only, the in- 
variance of concurrence under local unitaries, as well as 
the cyclic invariance of the trace operation, then make it 
possible to completely disregard the last rotation. The 
first rotation e~2Xi<^i ^ g-5X20-^ gjj^pjy translates into 
a relative phase X ■= Xi + X2- For an analysis of entan- 
glement and purity of a pure state under phase damping 
it is thus sufficient to study states of the form 
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